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AN EXPERIMENTAL INVESTIGATION OF TRANSONIC FLOW PAST 
TWO-DIMENSIONAL WEDGE AND CIRCULAR-ARC SECTIONS 
USING A MACH-ZEHNDER INTERFEROMETER 
By Arthur Earl Bryson, Jr. 

SUMMARY 


Interferometer measurements are given of the flow fields near two- 
dimensional wedge and circular-arc sections at zero angle of attack at 
high-subsonic and low-supersonic velocities. Both subsonic flow with 
local supersonic zone and supersonic flow with detached shock wave have 
been investigated. Pressure distributions and drag coefficients as 
fxmctions of Mach number have been obtained. The wedge data are com- 
pared with the theoretical' work on flow past wedge sections of Guderley 
,and Yoshihara, Vincenti and Wagoner, and Cole. 

It is shown that the local Mach number at any point on the surface 
of a finite three-dimensional body or an unswept 'two-dimensional body, 
moving through an infinite fluid, has a s-tationary val\ie at Mach manber 1 
and, in fact, remains nearly constant for a range of speeds below and 
above Mach number 1. On the basis of this concept and the experimental 
data, pressure, distributions and drag coefficients for the wedge and 
circular-arc sections are presented throughout the entire transonic 
reinge of velocities. 


INTRODUCTION 


Difficulties of Theory eind Experiment in Transonic 
Reinge of Velocities 


The difficulties inherent in studying transonic flow are well-known. 
Theoretical analysis is made difficult by the nonlinearity of the dif- 
ferential equations of compressible fluid motion. This nonlinearity 
leads to a change-o-ver in type of the differential equations from ellip- 
tic to hyperbolic when transition is made from subsonic to supersonic 
speeds. Since the essential feature of transonic flow is this mixed 
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subsonic- supersonic character, it is obvious that no linearization of 
the differential equations (at least in the physical plane) can 
adequately describe the flow. 

Wind-tunnel studies in the transonic range are made difficult by 
the large lateral extent of the perturbation flow field around bodies 
in this range. This means that models which are small compared with 
the test section must be used. Even then there is still a range of 
speeds from Just below = 1 to just above = 1 where the model 
and/or its support configuration are "choked,” that is, where local 
supersonic zones embedded in the subsonic field extend from the model 
to the tunnel walls, or, in the supersonic case, where embedded subsonic 
zones extend to the tunnel walls, or shock waves, reflected from the 
walls, impinge on the model. Some progress has been made recently in 
modifying wind-txmnel test sections so as to minimize these effects, 
but, on the whole, the majority of good test data in the range very 
close to M<jo = 1 has so far come from free-flight tests. Some good 
treinsonic data are available, however, from transonic-bump tests made 
in wind tunnels (reference l) . Using small models usually results in 
low Reynolds numbers so that difficulty is often experienced in extra- 
polating data to full-size Reynolds numbers; this seems to be particu- . 
larly true of the transonic speed range since the effects of boundary- 
layer and shock-wave interactions seem to be quite large there 
(references 2 and 3) • 

In this paper it is shown that in many instances tests need not be 
made in the region very close to Ma, = 1 since the flow in this range 
can be inferred from testing below and above this rsinge and using an 
interpolation based on the fact that the local Mach nianber at any point 
on the surface of unswept two-dimensional bodies and finite three- 
dimensional bodies has a stationary value at Mm =1. 


Existence of Potential Transonic Flows 

Guderley (reference U) has made a detailed investigation of the 
possibility of smooth transonic flows (i.e., subsonic flows with an 
embedded supersonic zone in which no shock waves appear) . He proposes 
that such smooth flows are exceptional, that they are discrete cases 
occ\irring for only particular body shapes at particular free- stream 
Mach numbers. Any perturbation of the shape with the Mach number held 
constant (or vice versa), Guderley claims, would result in a shock 
appearing in the flow. This bears an analogy to the well-known Busemann 
supersonic biplane which theoretically has ho shocks (and hence no drag) 
at a discrete value of free-stream Mach number and angle of attack 
(reference 5, p. 15^ • Guderley 's proposal is still controversial 
(e.g., see the paper of Sears who has made a critical survey of the work 
to date on the existence of transonic potential flows (reference 6)). 
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It is obvious that the potential flow must break down for a given, 
body shape at some Mach number less than .l. The argument whether this 
breakdown occurs precisely when a supersonic region first appears on 
the body or at a slightly higher Mach nijmber seems somewhat academic 
(although very interesting), since it is well-known experimentally that 
the drag-rise Mach number (i.e., the Mach number where noticeable shocks 
first appear) is very close to the critical Mach number (i.e., the Mach 
number at which sonic velocity first appears on the body) for most 
bodies without surface slope discontinuities. 

Kuo (reference 7) proposes that supersonic compression is unstable 
to disturbances; that is, a supersonic region on a body in subsonic 
flow must end in a shock with no compression occurring in the supersonic 
flow ahead of the shock. There seems to be ample experimental evidence 
to show that this is not strictly true since, for example, the compres- 
sion region of a X- shock is clearly supersonic. However, the X- shock 
configuration is believed to be a phenomenon associated with laminar- 
boundary- layer and shock-wave, interaction; with turbulent boundary layer 
(a condition more closely approaching nonviscous flow) hardly any' 
noticeable supersonic compression occurs before the shock ending the 
supersonic zone ( see reference 8) . 


Choice of Models 

Two-dimensional flow is much simpler to handle than axially sym- 
metric flow both in theoretical work and in interferometry. Hence it 
was decided to study two-dimensional flows despite the well-known dif- 
ficulties in approximating two-dimensional flow in a wind tunnel. 

Because of the considerations mentioned previously it was decided 
to test very small models which would be of such a shape that viscous 
influences would not materially affect the flow over them. This led to 
the choice of half airfoils ^ wedges and circular-arc sections followed 
by straight sections. These models have favorable pressure gradients on 
their surfaces over most of the transonic range so that boundary- layer 
separation, if it does occur, will only occur because of shock-wave 
influence. Furthermore, such separation will occur downstream of the 
part of the body being studied and hence will not affect the measurements. 
Certain viscous effects will still be evident, however, for instance, the 
effective rounding off of the shoulders and leading edges of the wedge 
models. 

Both theoretical advantages and practical need make the study of 
thin sections desirable. Consequently the semiwedge angles chosen 
1^ 1^ o o 

were ^ 10 (a 26.6 wedge was also used in order to make a 

comparison with some available theoretical work on a wedge of this 
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angle) . The circular-arc section chosen was essentially the front half 
of an 8 . 8 -percent- thick biconvex circular-arc airfoil, followed by a 
straight section. Models of sections much thinner than this, with the 
same chord lengths used, run into structural difficulties and also the 
ratio of boundary- layer thickness to model thickness becomes large 
enough to cause considerable deviation from nonviscous flow. 


Transonic-Flow Theory and Experiments 

The investigations of Von Ka^maji, Busemann, Guderley, Frankl, 
and many others have contributed significantly to methods of approach 
which can be used to study transonic flow (references 9 to l4) . The 
detailed numerical calculations for specific cases made by Maccoll and 
Codd, Emmons, Drougge, Drebinger, Guderley and Yoshihara, and Vincent! 
and Wagoner (references 15 to 21) have helped to dispel the idea of a 
"sonic barrier." Recently Cole at GALCIT has given an analysis of the 
flow past wedge sections at high- subsonic speeds (reference 22) . By 
combining the results of Guderley and Yoshihara 's, Vincent! and Wagoner's, 
and Cole's calculations, the flow past thin wedge sections can be given 
completely through the transonic range permitting a comparison with the 
present experiments. Some of the investigations mentioned above will 
be discussed in more detail further on in the present paper. 

Available experiments in the transonic range on thin wedge sections 
are surprisingly few. Pack (reference 23) describes some interferometric 
experiments on 10 and 20 semiangle wedges made at Braunschweig. His 
subsonic data appear to be good but the flow in the supersonic inter- 
ferograms appears to be very nonuniform and not very closely two- 
dimensional; only one supersonic Mach mmiber was tested where detached 
shocks occurred. His conclusion that the p/Poo distributions on the 
surface of the 20 semiangle wedge are very much the same for Mqo = O .803 
and Moo = 1.40 is interesting, but the statement that this agrees with 
the theoretical predictions of Maccoll and Codd is incorrect since they 
indicated that the p/pq distributions would be nearly the same. 

Griffith at Princeton has just recently published the results of 
some very carefully done experiments on flow past wedge sections of 
semiangles of 7°^ 10 °, 20 °, 30 °, ^ 5 °^ and 90 ° (and several other shapes) 
with detached shock waves (reference 24). These experiments were done 
in a shock tube and interferograms are presented of the flow fields. 

The experiments clearly show that the shape of the detached shock and 
its detachment distance from the sonic point on a wedge depend only on 
the body thickness and the Mach number (not the wedge angle) when the 
Mach number is well below the shock-attachment Mach number. This is in 
general agreement with Busemann 's considerations in his paper on detached 
shock waves ( reference 10 ) . 
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Liepmann, Ashkenas, and Cole (reference 8) made some careful pres- 
sure measiirements on the surfaces of 6- and 12-percent- thick "biconvex 
circular-arc airfoils at zero angle of attack at high-suhsonic speeds 
in connection with studies of shock-wave and boundary-layer interaction. 
Some of the results of their tests are combined here with corresponding 
low-supersonic test results from the present investigation to indicate 
the behavior of the pressure distribution on circular-arc airfoils at 
zero angle of attack through the entire transonic range. 


Acknowledgment 

While the present experimental investigation was in progress, a 
parallel theoretical study was being carried out by Mr. Walter G. 
Vincent! at the Ames Aeronautical Laboratory of the NACA. As pointed 
out in the introduction of reference 21, the two investigations were 
undertaken independently but the direct relationship of their subject 
matter made it desirable that they be reported similarly and it was 
therefore planned to present the results of both studies in a series of 
coordinated papers to be published by the NACA under a common general 
title. It has been decided, however, that the general title, "Transonic 
Flow past a Wedge Profile with Detached Bow Wave," does not completely 
describe the present report since high- subsonic- speed flows as well as 
low-supersonic-speed flows with detached bow wave were investigated and 
a circular-arc profile was tested as well as a wed^e profile. 

This work was conducted at the California Institute of Technology 
under the sponsorship and "with the financial assistance of the National 
Advisory Committee for Aeronautics. Grateful acknowledgment is made to 
Dr. Hans Wolfgang Liepmann who supervised the research, to the members 
of the GALCIT Transonic Research Group, to Mr. Walter G. Vincenti of 
the NACA, and to Dr. J. D. Cole of GALCIT for helpful discussions and 
assistance in carrying out the experimental work. Some of the results 
of these experiments have already been reported in reference 25. 


SYMBOLS 


a 

c 



sound velocity 

airfoil chord 

pressxire drag coefficient 

reduced drag coefficient 


\ (t/c) 5/3 ) 
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pressure coefficient 
reduced pressure coefficient 


/(7 * 
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Gladstone -Dale constant 

model span 

Mach number 

index of refraction 

pressure 

dynamic pressure 

airfoil thickness ratio 

horizontal component perturbation velocity (perturbation 
from a*) 

vertical component perturbation velocity 
Cartesian coordinates, origin at leading edge of profile 
reduced vertical distance f' + l)t/^^^^ 
angle of attack 

ratio of specific heats (l .4 for air) 
semiwedge angle 
wave length of monochromatic light used on interferometer 
reduced Mach number 
density 



Subscripts and superscripts: 

( ) conditions in free stream 

00 



reservoir conditions 
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( ) reservoir conditions behind a shock wave 

o 

{ ) conditions at sonic velocity 

Symbols used without subscripts indicate local conditions. 

APPARATUS AND METHODS 
Wind Tunnel 


The measurements were made in the GALCIT k- by 10- inch transonic 
wind tunnel. For a description of the tunnel and the flexible nozzle 
employed see reference 26. The tunnel can be run at both subsonic and 
low-supersonic velocities with continuous Mach number variation through 
use of the flexible nozzle and a variable second- throat nozzle down- 
stream of the test section. 


Models 

The models used were "half airfoils" followed by strai^t. sections. 
Foxir of the models were wedges (semiangles 4.53°, 7*56°, 10.00°, 
and 26.57°) followed by straight sections and the fifth was half of a 
biconvex circular-arc airfoil (8.80 percent thick) followed by a 
straight section (see fig. l) . The distance from the leading edge to 
the point where the straight section began was of the order of l/k inch 
for all five models. The models were made of tool steel and were very 
carefully machined and lapped so as to give exact cylindrical surfaces. 
Two pressure orifices on opposite sides of the airfoil were placed 
exactly the same distance from the leading edge to aid in setting the 
model to zero angle of attack by balancing these pressures on an alcohol 
U-tube. Because of the very short chord lengths vernier-protractor 
measiirements of the opening angles of the leading edge were of doubtful 
accuracy, so the angles were measured by letting the leading edge split 
a beam of parallel light and measuring the position of the reflected 
spots on a wall in back of the model. In this manner the angles could 
be measiired to ±0.03°. 


Interferometer 

The interferometer used in this investigation is described in 
references 27 and 28. One of the main features of this interferometer 
is that both li^t beams are passed through the test section, one over 
the model and the other ahead of the model in the uniform flow field, 
that is, where the velocity is nearly the free- stream velocity. The 
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advantages of this are: (l) The fringe shifts are in relation to the 

free-stream density and (2) the effects of the side-wall hoxmdary layers 
are approximately canceled out since both beams traverse nearly the same 
boundary layer at each side window. This leads to improved accuracy 
when the interferograms are evalviated on the basis of the absolute value 
of the fringe shift from no-flow conditions. For these tests finite- 
fringe interferograms were used and another method of evaluation was 
devised which is much simpler and more accvirate than the above-mentioned 
technique. Infinite-fringe interferograms, while they give the constant- 
density contours immediately, are less accurate than the superimposed 
finite-fringe interferograms because any optical inaccuracies in the 
system cause the contour fringes to be distorted. These inaccuracies 
are calibrated out in the superimposed finite-fringe interferograms. 

Also there are times when one does not know whether the density incre- 
ment between contours of an infinite-fringe interferogram is positive 
or negative; this trouble does not arise with the finite-fringe inter- 
ferograms. A typical finite-fringe interferogram is shown in figure 2. 


Method of Evaluation of Interferograms 

The method of evaluation used here depends on two techniques: 

(a) Photographic superposition of disturbed and undisturbed interferograms 
and (b) fringe identification by a pressure measured on the model. 

Direct photographic superposition of a "no-flow" finite-fringe 
interferogram on a "with-flow" finite-fringe interferogram gives rise 
to dashed shadowy lines (the dashes being where the dark fringes of 
one picture cross the light fringes of the other) ; see figure 3 for an 
example of this type of picture. These shadowy lines can easily be 
shown to be lines of constant density for two-dimensional flow and are 
the same contoxirs as would be obtained on an infinite -fringe inter- 
ferogram made with perfect optical surfaces. The increment in density 
between these shadowy lines is a constant dependent only on the span of 
the model and the wave length of the monochromatic li^t being used. 

This is easily shown since the difference in optical path lengths of the 
light rays between two adjacent constant-density contours must be 1 wave 
length of the light being used. For two-dimensional flow the difference 
in optical path length will simply be I An, where I is the span of • 
the model and An is the difference in index of refraction between the 
two light paths. Thus 


I /Sn = X (1) 

But the relation between index of refraction and density in a gas is 
given by 


n - 1 = kp 


( 2 ) 
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where k is the Gladstone-Dale constant (a function of the light 
frequency and type of gas) . Therefore 

( 3 ) 

where Ap is the difference in density between two adjacent constant- 
density contovirs. For these experiments 

k = 5^61 A (mercury green line) 

k = 0.1162 cubic foot per slug 

I = 3*50 inches 

so 

^ = 0.0250 per fringe shift 
^o 

where p = 0.00211 slugs per cubic foot was the usual tunnel stagnation 
density.® 

The advantage of photographic superposition is not only in time 
saved but also in increased acci^acy of evaluation. Any slight changes 
in fringe spacing or fringe orientation with respect to the "no-flow" 
interferogram which occur before the "with-flow” interferogram is taken 
can be almost exactly canceled out by causing the two superimposed 
interferograms to coincide exactly in a region where it is known that 
the flow was imiform, since in such regions there should be no isopycnic 
contoiars . This is particularly easy to do for supersonic flow if a 
portion of the flow field ahead of the nose shock wave is included in 
the interferogram. For subsonic flow care must be taken to include 
enough of the flow field ahead of the model in the interferogram to 
have some of the nearly undisturbed flow field for comparison; this was 
quite simple to do for the small, thin models used in these tests. 

The actual superposition technique used here was first to a 

print times enlarged) of the "with-flow" interferogram. This print 

was then placed under the enlarger and the "no- flow" interferogram 
negative was put into the enlarger. By changing the enlargement scale 
and moving the "with-flow" interferogram under the enlarger the fringes 
were made to coincide exactly in the regions of uniform flow. The 
constant-density contours could then be drawn in on the print. Alter- 
natively the first print could be made on transparent paper (Ansco 
Reprolith Ortho was used) and when the superposition was accomplished a 
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piece of photosensitive paper was slipped under the transparent print 
and a print of the two interferograms was obtained. This was the 
technique used for figure 3 - 

• In Order to identify the density values with the fringes a pres- 
sure tap was placed on the model in a region where the pressiire gradi- 
ent was expected to he large. From the pressiare reading the density 
at the pressure tap was calculated using reservoir fluid properties ' 
(taking into account entropy changes through shock waves). The pressure 
tap will always lie between two fringe contours or on a contour, so that, 
by knowing the density increment between fringe contoiars, the values of 
the density on the adjacent contours can be obtained by interpolation. 

The whole interferogram is determined once the density is known on one 
contour (except for the shifts through shock waves). 


Interferometer Sensitivity 

It is interesting to note that the interferometric method has its 
greatest sensitivity in the transonic range. As pointed out, previously, 
the density increment between two adjacent contour fringes is a constant 



Now in any part of the flow field where the stagnation density is con- 
stant along a streamline. 


so 



ih) 




M dM 


( 5 ) 


Hence the increment in Mach number between adjacent contour fringes is 
given approximately by 



kZpo 


AM.~ 


M 


( 6 ) 
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since 

Ap _ \ 

p ~ kZp 
o o 


This function has a minimum at M = y- ^ which is M = , 0 . 9 l 4 for air 

A graph of this function is shown in figure 4 . Note that 
the Mach number increment per fringe for these tests was always closely 
equal to 0.05 . 


Similarly, the expression for the increment in pressure coefficient 
between adjacent contour fringes is approximately 


AS 






7 





X 

kZPo 


( 7 ) 


For values of M close to this expression has a minimum at 

=\~2~^ which is I.832 for air. 


Side-Wall Boundary-Layer Effect on Approximating 
Two-Dimensional Flow 

A close approximation to two-dimensional flow over the whole span 
of the model was required since the interferometer integrates the value 
of the density from wall to wall. In a nonviscous fluid, letting the 
model extend from wall to wall would theoretically give two-dimensional 
flow over the whole span. If the model did not span the whole tunnel, 
the flow would correspond to that past a model of infinite span with 
periodic gaps in it where the gaps were equal to twice the distance from 
the edge of the' model to the wall.* The effect of the side-wall boundary 
layers, for a model that does not span the tunnel, is rou^ly to decrease 
the size of this gap. Approximately, the gap size would be decreased by 
twice the displacement thickness of the wall boundary layer. By making 
the gap between the edge of the, model and the wall approximately equal 
to the wall-boundary- layer displacement thickness, one might hope to 
approximate closely two-dimensional 'flow over the span. This phenomenon 
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is, of course, very much more complicated than this, particularly in 
the supersonic case where the shock waves interact with the wall bound- 
ary layer. However, by taking circular-cylinder and wedge models and 
varying the gap size in increments of 1/16 inch, it was found that the 
detached bow wave became closely two-dimensional when the gap size was 
1/4 inch (i.e., there was no blur ahead of or behind the shock pictxires) 
which is almost exactly the boundary- layer displacement thickness when 
measured without a model in the test section. When the gap was 3/16 inch 
the shuck was blurred ahead of the main shock and when the gap was 
5/16 inch it was blurred behind the main shock. These tests were fur- 
ther substantiated by §ome schlieren pictures, which Mr. Walter G. 
Vincenti of the NACA Ames Aeronautical Laboratory kindly made available, 
showing a view looking down on a wedge model so that the leading edge 
of the detached shock appeared as a line; by varying the model span a 
discrete value of the span was found where this line was almost exactly 
parallel to the leading edge of the model, while for just slight vari- 
ations from this gap size the shock was curved forward or backward. 

Figure 2 shows a finite-fringe interferogram of the circular-arc section 
with a detached shock where the definition of the shock wave was 
unusually sharp. This is strong, but of course not conclusive, evidence 
that the flow was closely two-dimensional over most of the span. Fur- 
ther evidence that the flow differed from two-dimensional flow only 
slightly is given in the next section. 


Side -Wall Boundary-Layer Effect on Interferogram Evaluations 

A result of the method of interferogram evaluation described above 
is that the effect of the side-wall boundary layer is approximately 
canceled out, since the over-all fringe shift from no-flow conditions 
is unimportant^ only the relative fringe shifts from a point of known 
density being used. This is strictly true only if the integrated side-- 
wall boundary- layer density, defined by 

6 

P dy (8) 

where y is the direction perpendicular to the tunnel wall and y = 0 
is the wall, is the same over the entire field of view of the interfer- 
ometer. Obviously, this can never be exactly true since the pressure 
field caused by the model, the boiondary- layer growth, and the shock- wave 
and boundary- layer interaction all tend to change this value. An indi- 
cation that all these effects might be small was obtained from the model 
tests where pressures were measured at two points on the model in the 
center of the span, where the flow is closely two-dimensional; the den- 
sity increment between these two points on the model was compared with 
the density increment given by the interferogram. The standard deviation 
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from zero of the difference between these two increments over the whole 
range of test Mach numbers was about 1 percent of the stagnation den- 
sity. Also, the values of pressure-drag coefficient obtained inter- 
ferometrically for the attached-shock-wave- cases checked the oblique- 
shock theory very closely, and it is well-known that the oblique-shock 
theory checks experiment quite well. 


Determination of Free-Stream Mach Nvunber 

An interesting result of the method of evaluation Just described 
is that the free-stream Mach number in subsonic flow can be determined 
from the interferogram and the measured pressure on the model, provided 
a large enough field of view ahead of the model is obtained in the 
interferogram. This can be done by noticing that a certain number of 
compression contovirs appear around the leading edge and then expansion 
contours follow these toward the back part of the airfoil; the center 
fringe corresponding to free-stream density can then be traced out into 
the flow field (see, e.g., figs. 9(a) to 9(d) of the 10° wedge in sub- 
sonic flow). The exact value of the density can be determined on this- 
fringe as described previously and, hence, knowing the stagnation den- 
sity in the settling chamber, the effective free-stream Mach number can 
be determined from the isentropic-flow relations. It is believed that 
this effective Mach number is a good approximation to the free-flight 
free-stream Mach number and would give the same flow as that measured in 
the wind tvinnel for the very small models used in these tests. 

This method is more accurate at high- subsonic speeds than at low 
speeds since more contour lines are obtained on the airfoil at the higher 
speeds (see above discussion). The estimated accuracy in determining 
free-stream Mach number in this way was ±0.01 for the range of subsonic 
Mach numbers tested. 

The free-stream Mach ninnbers for the siq)ersonlc tests were obtained 
by calibrating the, flexible -nozzle Jack settings against Mach number 
with a static-pressure probe in the center of the tunnel. The probe was 
traversed upstream and downstream in the region where the models were to 
be tested and an average Mach number was obtained there. The standard 
deviations from this average value were of the order of ±0.005 in Mach 
number for the range of supersonic Mach mambers tested. 


Wind-Tunnel Choking 

In all the subsonic testing the embedded supersonic zone was not 
allowed to touch the upper or lower walls. In one or two of the low- 
siqjersonic tests there was a question whether the embedded subsonic zone 
touched the ceiling or not. In case it did, it is well-known that in 
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such cases the detached shock changes its curvature near the ceiling so 
as to come in nearly normal to the walls. Since the models were so very 
small (1/16 in. thick .compared with the 10-in. height of the tunnel), it 
is believed that the effect of this on the pressure distribution was 
negligible. 


Reynolds Number 

The value of the Reynolds number for all of these tests was approx- 
imately 60,000 based op the chord of the model. The boundary layer on 
the models was laminar and no effort was made to trip the boundary layer 
to make it turbulent. The compression region in the shocks shown in the 
high- subsonic-flow interferograms is believed to be associated with the 
laminar boundary layer, as mentioned previously. 


THEORETICAL WORK ON TRANSONIC FLOW 
Relaxation Calculations 


In 1946 Maccoll presented a paper at the Sixth International 
Congress for Applied Mechanics in which he described a relaxation cal- 
culation of the compressible flow past a 20° semiangle wedge followed 
by a straight section at Mach ntimbers of O.7 and 1.5. The flow field 
in both cases contained both subsonic and supersonic velocities. His 
main assumptions were; (l) Sonic velocity occurs at the shoulder and 
(2) the streamlines of the flow are perpendicular to the sonic line 
(i.e., the line where sonic velocity occurs in the flow). The first 
assijmption can be shown to be correct (see reference 22) so that indeed 
it is not an assumption. The second assumption, as Maccoll realized, 
was only approximately correct for M^o = 1.5 and certainly quite incor- 
rect far away from the wedge at M„ = 0.? (since the asstunption leads 
to an infinite supersonic region above the wedge) . In effect his solu- 
tion at Moo = 0*7 was "choked" in the sense that the back part of the 
body could have no influence on the front part. It is well-known that 
for bodies , at high-subsonic speeds a finite, closed supersonic region 
occurs in the flow, so that the sonic line makes all angles possible 
with the streamlines, including 0 °. The method of solution used was to 
assume positions of the shock wave and sonic line, calculate the resi- 
dues in the relaxation net using the isentropic-flow equations (an 
approximation since flow behind a curved shock is not isentropic) , and 
then readjust the shock-wave and sqnic-line location, calculate again, 
and so forth, iterating until the solution closely repeated itself. 
Maccoll found that the p/pq* distribution on the wedge surface at 

= 1.5 was nearly identical with the p/p^ distribution at M„ = 0.7. 
This led him to propose that the pressure in the transonic region, on 
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todies with distinct corners, varied as the stagnation pressure and he 
presented a drag curve through = 1 for the 20° semiangle wedge cal- 
culated on this basis. 

Drougge in 19^*8, following Maccoll, calculated the flow past a 
finite cone of semiangle with' detached shock wave at . = I.80 

and = 2.15, using the same assumptions as Maccoll (reference 18) . 

He also made experiments on this cone and found the agreement with his 
theory rather good. He made se'veral tests at lower supersonic Mach 
numbers also and found that the p/Pq* distribution on the cone surface 
did remain nearly constant except as the Mach number became close to 
the attachment Mach number. 

Drebinger in 1950 showed how to calculate, by relaxation techniques, 
the flow past finite cones and wedges with detached shocks, eliminating 
the isentropic-flow assmption and the asstunption on the streamlines • 
being perpendicular to the sonic line (reference 19) . He calculated a 
specific example - a 26.6° semiangle wedge at M^ = 1.44o - and checked 
the calculated shock-wave shape and position experimentally. His cal- 
culations showed that, even for the detached- shock case, the streamlines 
differed from being perpendicular to the sonic lines by angles as large 
as 30 . His calculation was checked in detail experimentally in these 
tests and agreement was found to be excellent. 


Transonic Perturbation Theory 

By assuming that the velocity component parallel to the free stream 
direction differs only by a small quantity u from a*, the critical 
velocity, and keeping only the highest-order terms in the differential 
equation, the equations of two-dimensional irrotational fluid motion are 
reduced to 


(7 + 1) 


u 8u 8v, 
a* ^ 




= 0 


8u ^v _ ^ 

J 


(9) 


It was from these equations that Von Karman and Guderley independently 
arrived at the transonic similarity laws (references 9 and k) , For two- 
dimensional steady flow past sections whose shape fiinctions are the 
same, these laws imply that . - 
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M^ - 1 

- f, 

< 


-1 

Wi 


Wi' 

(7 + l)t/c 


(7 + l)t/c 

. -> 


where M is the local Mach number on the surface of the section, 
similarity in pressure and drag coefficients is then 


( 10 ) 


The 


(7 + 1)^'^% 


(t/c) 
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(7 + 


(t/c) 


ITT 


= h 


r 

- 1 


J 

1 1 

+ 

rt- 

0 

1 1 


1 — 1 

1 

(^J 

.L 

(7 + l)t/c 

p/3 

1- _ 
V 



( 11 ) 


( 12 ) 


These quantities will be called reduced local Mach number, reduced free- 
stream Mach number, reduced pressure coefficient, and reduced drag coef- 
ficient, respectively, using symbols 5, | , C , and 

00 p 1 ) 

By interchanging dependent and independent variables in the 
perturbation equations the problem becomes linear 


“ ^ = o' 

ov ou 

^ ^ = 0 

' ov c3u 


V 


(13) 


where 


= (r -M) ^ 


V = (7 + 1) 


a^ 


and by eliminating x by differentiation, the Tricomi equation is 
obtained 
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The main difficulties with this holograph (u,v) plane are: (a) The 

mapping^ of physical boundaries into the holograph plane is in general 
not known until the solution to the problem is known so that it is not 
known where to apply the boundary conditions in the holograph plane and 
(b) the mapping is often multivalued, complicating the solution. Two 
interesting cases are known where these difficulties are avoided. They 
are: (a) The free jet, studied by Tschaplygin in 1905, and (b) the 

finite wedge, studied recently by Guderley and Yoshihara, Vincent! and 
Wagoner, and Cole. These latter studies came to the author's attention 
after the present experimental study of the finite wedge in trdnsonic 
flow had begun and served to make the study more interesting since the 
data could then be compared with the theoretical results. 


Theoretical Studies of Transonic Flow past 
Thin Wedge Sections 

Guderley was apparently the first to formulate the problem of the 
thin finite wedge in the holograph; he and Yoshihara fovmd an approximate 
solution to the problem of the flow past a thin double-wedge profile at 
zero angle of attack at Mach number 1 using the transonic perturbation 
equations (reference 20) . 

Vincent! and Wagoner considered the thin double-wedge profile at 
zero angle of attack for low- supersonic reduced Mach nvmbers where the 
shock wave is detached (reference 21) . Their solutions were effected 
by relaxation calculations in the holograph plane. Here the bow shock 
wave and the sonic line are fixed boundaries (their positions are not 
known originally in the physical plane) and the boundary condition on 
the shock is the slope of the streamlines (or y = Constant lines). This 
boundary condition was first shown by Busemann, , who aptly called the 
configuration a "hedge hog." 

Cole (reference 22) has recently given a simple approximate 
analytical solution to the flow past a thin wedge at zero angle of 
attack followed by a strai^t section at high- subsonic speeds (M^j, < l) . 
His solution satisfies the Tricomi equation and the boundary conditions 
on the wedge and at infinity but not the boundary conditions on the 
sonic line. Effectively his solution gives a finite vertical sonic line 
from the shoulder which is also a limiting line. Cole has indicated 
that this solution is the singular part of the solution in the hodograph 
and as such is most likely the main part of the solution. It is 
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interesting to note that the drag-curve slope and cvirvatxire at Mm = 1 
obtained from Cole's solution agree exactly with the values obtained 
from the simple physical considerations of the next section. Also, the 
pressure distribution on the wedge at Moo = 1 agrees within 1 or 2 per- 
cent with that obtained by Guderley and Yoshihara. 

Since the back half of a double-wedge profile has only a very weak 
influence on the pressiare distribution on the front half for > 1 
(only through the "last Mach wave" from the shoulder point to the sonic 
point on the detached shock) , it is reasonable to take the solution of 
the double wedge at M^, > 1 and use the front-half solutions in con- 
nection with Cole ' s results for M^^ < 1 for the wedge followed by a 
straight section and thus have a solution for the latter semi-infinite 
body completely throu^ the transonic range. By using linearized sub- 
sonic theory and the shock-expansion supersonic theory, the zero-angle- 
of-attack flow is obtained for all possible values of 

Tsien and Baron (reference 29) have shown that the shock-expansion 
theory can be expressed in the transonic similarity form for thin bodies 
in piore supersonic flow near M^^ =1. 

Von KaWan (reference 9) has Indicated also how linearized subsonic- 
and supersonic-flow resvilts may be written in the transonic similarity 
form since, from the Prandtl-Glauert similarity, in linearized subsonic 
theory, 

1) (^ 5 ) 


and, in linearized supersonic theory. 




and from the expressions for reduced pressure coefficient and Mach nuro- 

(r + 

her, multiplying both sides by — — ^ — o 7 ~ > "these equations may also be 

(t/c)^/" 


written as: 



(15a) 
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Moo^ - 1 



so equations ( 15) and ( 16) may be written in transonic form 



The subsonic pressure-distribution and drag-coefficient curves 
have been calculated here from Cole's analytical expressions and, com- 
bined with the results of Guderley and Yoshihara, Vincent! and Wagoner, 
and Tsien and Baron, the curves for reduced pressvire and Mach number 
distribution and reduced drag coefficient^ are given in figures 5 to 8 
for the finite wedge followed by a straight section. 


^he reduced drag coefficient given in figure 8 is that for the 
half wedge and is equal to 
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It can te shown that Cole's solution for large negative values 
of goes over exactly into the linearized subsonic solution (see 

appendix A) . The reduced-pressure-coefficient curve for = -2.02 
in figure 7 is so nearly identical for both solutions that they cannot 
be told apart (except that Cole's solution goes to C^ at x/c = 1 

while the linearized solution goes to -«) . This is to be expected 
since the transonic perturbation equations are not restricted to tran- 
sonic flow but apply equally well to completely subsonic and completely 
supersonic flow. 2 The transonic equation can be written in the form 




3x2 


(19) 




where <p is the perturbation potential such that u = U + 

Thus it is clear that for completely subsonic or completely supersonic 
flows the term on the right is negligibly small but becomes of paramount 
importance in transonic flow. 


CHARACTERISTIC FEATURES OF TRANSONIC FLON PAST 
WEDGE AND CIRCULAR-ARC SECTIONS 
Characteristic Free-Stream Mach Numbers 


Critical Mach number .- The Mach number at which sonic velocity 
first appears on the wedge is Moo = 0 (within the inviscid theory) 
since subsonic flow can not t\im a sharp comer. Because of the fact 
that the botmdary layer rounds off the comer, and perhaps also because 
of the spatial resolution limitations of the interferometric method, 
sonic velocity was not found there experimentally until approximately 
= -O. 8 O for the wedges. 

The critical Mach number for a half circular-arc airfoil followed 
by a straight section can be obtained approximately from linearized 
subsonic theory. This theory gives the surface pressure distribution as 



^his was pointed out to the author by Dr. Milton Van Dyke of the 
NACA Ames Aeronautical Laboratory 
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which yields ' 


■min 


^ -1.626(t/c) 


( 21 ) 


M 


at x/c = 0.783 (see appendix B) . This equation can also be written 

(r + 1)^/3 

in transonic similarity form by multiplying both sides by — — ^ - ■ 

(t/c)^/3 

( as shown in the previous section) : 


■ min 


■1.626 


(22) 


Now, within the transonic perturbation theory 


Cp ■= -2(S - U 


( 23 ) 


Hence 


^p “ ^^00 

■^cr 


( 24 ) 


Equating C_ to one obtains the critical reduced Mach number 

■^min ■^cr 


loo = -0.871 

cr 


For the thickness ratio t/c = O.O88 used in these tests, this predicts 
a critical Mach nimiber of 0.834 at x'/c = 0.783* Experimentally the 
critical Mach number was found to be 0.825 and occurred somewhere 
between x/c = 0-75 and 0 - 95 _(the pressure distribution was very flat 
in this range) . It is interesting to note that the experimental M 

®“cr 

was higher for the wedges than for the circular-arc profile of the same 
thickness ratio. This was probably due to a combination of three 
effects: (l) The boundary layer for the same Reynolds numbers used 

here was fairly thick in comparison with the dimensions of the model 
and thus it "rounded off" the shoulder more than would be the case at 
higher Reynolds niambers. ( 2 ) The height of the supersonic zone, even 
for an ideal nonviscous flow past thin wedges, appears to be quite 
small until the free-stream Mach number is quite close to 1 . This is 
apparent from Cole's theory and also from the argument in reference 25 
that the height of shocks in the supersonic zone must be of the form 
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(r 1) ^ 


nV3 



(3) The spatial resolution of the interferometric method may not have 
heen sufficient to detect very small supersonic zones near the shoulder. 
There is also a large refraction error near the shoulder due to the 
high density gradients which tends to obscure details of the flow there. 


Shock-attachment Mach number .- The shock-attachment Mach number 
depends only on the opening angle of the profile at the leading, edge n 
and can be predicted quite precisely by oblique-shock theory. If 0 is 
the semiopening'' angle then it can be shown that approximately, for thin 
profiles, 




(7 + i)e 


Tfl ~ 



(25) 


(see appendix C) . If t/c is the thickness ratio of the circular-arc 
section, 0 2(t/c). Hence for the circular-arc profile 



(26) 


Mach number at which sonic velocity appears behind an oblique shock . 
The Mach number at which sonic velocity appears behind an oblique shock 
Moog is' just slightly higher than Ma,^ and again is a function only of 

the opening angle. These values can also be found quite precisely from 
oblique-shock theory and approximately in similarity form can be given 
as 


- 1 


(7 +. I)e 


2/3 


= 2 


1/3 


for the wedge (see appendix C) and 


I '= 2 

Soog 


(27) 


(28) 


for the circular-arc section 
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Characteristic Values of Local Mach Number 

Mach number at leading edge . - The Mach number at the leading edge 
is zero (a stagnation point) for all free-stream Mach numbers less than 
the attachment Mach nvnnber. 

Mach number at shoulder of vedge . - The Mach number at the shoulder 
of the wedge just before the turn is always 1. This is easily seen in 
the case of flow with detached shock since the only characteristic 
distance of the finite wedge is the distance from the leading edge to 
the shoulder which must determine the shock-detachment distemce^ and, if 
the sonic point occurred ahead of the shoulder, the shoulder could not 
influence the shock position. Subsonic flow cannot turn a shaip) corner 
so the flow must therefore reach Mach nvnnber 1 right at the comer. In 
the case of subsonic free-stream flow the argument is not so simple 
( see reference 22) . 

At the shoulder the flow around the comer is locally a centered 
Prandtl-Meyer fan starting from M = 1. The Mach number Just behind 
the comer is thus determined only by the wedge angle and is independent 
of the free-stream Mach nvnnber. Behind this point the flow will recom- 
press to the free-stream Mach nvnnber throu^ a shock or series of shocks, 
for free-stream Mach numbers less than the attachment Mach number. The 
expression for Mach nvnnber Mpjj behind an expansion from M = 1 

through an angle 6 is 

® 1 - 1 (29) 

Expanding the right-hand side in powers of - 1 (assvimed small) 

the first nonzero term yields 


e 




which is in transonic similarity form so 


(30) 



(31) 
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PRESSUBE DISTR3BOTI0N ON BODIES MOVING THROUGH AN 
INFINITE FLUID AT SPEEDS NEAR MACH NUMBER 1 
Stationary Value of Local Mach Numher at Free-Stream 

Mach Numher 1 


Diiring the course of these investigations it was found that for 
the wedge and circular-arc sections the local Mach numher distributions 
on these sections at very high subsonic speeds |ahove hut below 

choking Mach niomher) and at very low supersonic speeds (where the 
detached shock wave was a chord length or so ahead of the section) were 
nearly identical. In trying to understand why this should he so, the 
following explanation was derived: (l) At low- supersonic speeds the how 

shock wave is detached a great distance ahead of the profile and a sub- 
sonic flow region is embedded in the flow field between the shock and 
the sonic line. The part of the shock directly ahead of the profile is 
nearly normal over quite a distance (of course, the slope of the shock 
asymptotically tends to the slope of the Mach wave of the free- stream 
flow at large distances lateral to the flow direction) . Nagamatsu 
(reference 30 ) has previously indicated this and points out that. the 
flow past the profile should be closely approximated by assuming the 
profile, is in a high-speed subsonic flow where the velocity distribution 
at infinity is slightly nonuniform, the minimxaa velocity being directly 
ahead of the profile and equal to the velocity behind the normal shock 
and then increasing in both lateral directions. (2) Now the normal 
shock near Mach nimiber 1 is nearly symmetrical in the sense that the 
Mach number behind the shock is just as much below 1 as the Mach number 
ahead is above 1. This follows from the normal-shock relation: 



(32) 


where M]^ is the Mach n\jmber ahead of the shock and M 2 is the Mach 
number behind the shock; so near M^ = 1, 


2 2 
1 - Mq M^*^ 


( 33 ) 


or 


1 - Mg « M^ - 1 


m 
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Therefore if = 1 + e where € is small, the flow past the profile 
is nearly the same as the flow past the profile at = 1 - e since 

the Mach number behind the central part of the detached shock wave is 
almost exactly 1 - e. It follows therefore that the local Mach number 
distribution on the profile surface must have a stationary value at 
Moo = 1 and furthermore vary only slowly in the neighborhood of Moo = 1. 
Mathematically this means 


dM 

dMoo 


= 0 


M =0 


( 35 ) 


It shoiild be noticed that this argument is based on two assump- 
tions:^ (1) The detached bow wave moves very far ahead of the profile 
as the flight Mach number decreases toward 1 . ( 2 ) The radius of c\arva- 

tvire of the detached bow wave at points directly ahead of the profile 
becomes extremely large as the flight Mach number decreases toward 1. 

Examining these assumptions it would seem that the same reasoning 
should apply to any finite three-dimensional body in an infinite fluid 
traveling at speeds neax Mach number 1, except that now two radii of 
ciorvatxire at points on the detached bow wave ahead of the body must be 
assumed to become large as the flight Mach number decreases toward 1. 

The detached bow wave is so far away from the body at speeds Just 
slightly above Mach number 1 that the body appears as only a very small 
object in relation to the radii of cxirvature of the bow wave and, hence, 
it would appear as thou^ the shape and attitude of the body could have 
no appreciable effect in changing the argument presented above. 

The reasoning should also apply to an infinite yawed cylinder 
(whose cross section may be finite or, if the angle of attack is 0, may 
extend infinitely far downstream) provided that the Mach number con- 
sidered is the component of the Mach number normal to the generators of 
the cylinder. 

These arguments are for steady-fli^t speeds. Large accelerations 
throu^ sonic flight speed could conceivably modify the phenomenon. 

Thus it is difficult to Judge whether or not the available flight-test 
data confirm the concept since nearly all such data come from missile 
tests that Involved large accelerations (or decelerations) through sonic 
flight speeds. The transonic -bump tests of Weaver on sweptback wings 
(reference l) would seem to support these conclusions since they show 

3 , — 

It is believed that these are not actually assumptions but are 

capable of demonstration if one assumes a smooth variation of drag 
through Mach number 1. 


26 


NACA TN 2560 


drag-coefficient maximums very near Mach number 1, a necessary con- 
sequence of the concept for finite three-dimensional bodies and finite, 
unswept, two-dimensional bodies as will now be shown. 


Slope of Pressure- and Drag -Coefficient C\irves at = 1 

Equation (35) enables one to calculate the slope of the pressure- 
and drag-coefficient curves at Mach nianber 1 as follows: 



(36) 


7M 


/l . 2_^ 




- 1 


for isentropic flow so 


dC„ 

P 


2Cp| 


M«=l 


M„=l 


7 + 1 


7 + 1 


using 


(37) 


(38) 


dM 


= 0 




How for a two-dimensional body the pressure-drag coefficient (based on 
the chord) is given by the conto\ir integral 




t . S) ds 


(39) 


where 

i unit vector in stream direction 

n unit vector nozmial to profile pointing outward 

ds element of length along profile contoiir 
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so if the angle of attack is constant and Moo is changing 


dCjj 

dMoo 


M„o=l 



2Cr, 


P 



rrr-j 


^ -A. 

i • n ds 




But 



n ds = 0 for a closed contour, so 


<3Moo 


M =1 
00 


7 + 1 


Cd 


Moo=l 


ihl) 


For the front part of a profile (defined as 
maximiun thickness) the usual definition of a drag 


that part ahead of the 
coefficient is 



(42) 


vhere 


means the counterclockwise line integral from the point of 


a 

maximum thickness on the upper surface to the point of maximiim thickness 
on the lower surface; thus 


so 


dMoo 


Moo=l 



2C 


^ M =1) 


7 + 1 


i 


n ds 


where 


dC 


% 


dK 


Moo=l 


4 t 

TTT c - 


7 + 1 


"Dp 


M =1 
00 


( 43 ) 


t maximum thickness of profile 

a angle of attack of profile 
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Similarly the drag coefficient for the rear part is 





so 


dC 

dM 


% 




' k t 

— cos a - 

7 + 1 c 


7 + 1 




( 45 ) 


For the tests on wedge and circular- arc sections followed by 
straight sections the concept of drag coefficient of the front part of 
the section will often be used. 

For bodies of revolution (which includes spheres, cone- cylinders, 
etc.) the pressure-drag coefficient (based on maximxan cross-sectional 
area) at zero angle of attack is 


where 

R maximum radius of body 

I length of body 

X distance from nose along axis 

Therefore 


dCD 


-2Ct 


M„=l 


M =1 

00 


7 + 1 


as before in the two-dimensional case. However, for front and back drag 
coefficients 



( 47 ) 
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so 


and similarly 


dC 


dMco 


M„=l 


7 + 1 7 + 1 


Cpp 


M =1 
00 


(l*8a) 



<3M„ 


Mco=l 


4 

7 + 1 


. 2 

7 + 1 



(48b) 


and these differ from two-dimensional valiies obtained above in equa- 
tions ( 43 ) and ( 45 ) by not involving the fineness ratio of the body 
(this is of course due to the different reference areas for drag 
coefficients) . 


For the general finite 
coefficient is given by 


three-dimensional body the pressure-drag 


JJ (*^ 9 ) 

S 

where A is some reference area of the body and S is the surface of 
the body. It follows as it did previously that 


dCjj 

<3Mco 


^ 00=1 


7 + 1 D 


M =1 


Slope of Drag-Coefficient Curve at = 1 in Transonic 
Similarity Parameters for Two-Dimensional Flows 
Within the transonic approximation 

= -2(1 - ^„) (50) 

so 



-2 



( 51 ) 
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Now 


dM 

dMo 


= 0 implies that 


M =1 
*00 




L=0 


= 0; hence 


dCj 

dF 


= 2 , 


^ =0 


Now 


= "t^ i ) 


(52) 


(53) 


SO 





Similarly it is easy to show that 


dC 







and 


m 


(55) 





(56) 


other Data Showing Slow Variation of Local Mach 

Number Near M = 1 

00 

As mentioned previously, Maccoll in 19^6 had already proposed the 
slow variation of local Mach number near =1 on "bodies having 
distinct comers," It appears that this latter restriction is not neces- 
sary. Maccoll *s proposal was based on rather slim evidence and it is 
believed that here, on the basis of the argument presented concerning 
the normal shock, the principle is explained more convincingly. Also 
the experimental evidence given here and by Drougge (reference l8) , 
Bleakney and Griffith (personal communication), Weaver (reference l) , 
and by some NACA reports tends to bear out the conclusions of slow varia- 
tion of local Mach number on bodies near Moo = 1. 
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This fact is sometimes slightly obscured in the NACA reports 
because pressure coefficient was plotted instead of p/po or local 
Mach number. However; constant Mach number lines were sometimes drawn 
in these plots and there the evidence shows up strongly (see, e.g., 
reference 31> figs^ 7 to 11, pp. 36 and 37)*, The relative constancy of 
local Mach number distribution near M„ = 1 for airfoils at an angle 
of attack is also shown clearly in figures 8 , 9 , and 10 of reference 32 . 


ON COMPARING THEORY AND EXPERIMENT 


In references 21 and 25 discussions were presented on the philos- 
ophy of comparing experiments with approximate thories, and these dis- 
cussions will not be repeated here, except to mention that in some of 
the theoretical curves presented here the values have been shown with a 
certain spread which results from using a pressure coefficient equal 


to -2 


u - u 

IF 


00 


or -2 


u^ 

U 


(the former value is the one that fits 


into transonic similarity theory; the latter value is the one more com- 
monly used in perturbation analysis) . 


In connection wxth the idea presented in reference 25 of extra- 
polating experimental data to zero thickness in order to compare with 
results from transonic perturbation analyses, it is interesting to note 
that the characteristic Mach nvnnbers mentioned in the section "Charac- 
teristic Features of Transonic Flow past Wedge and Circular-Arc Sections" 
can be presented in powers of the thickness of the wedge (or equivalently 
in powers of the wedge angle), the first term of which gives the tran- 
sonic similarity expression; two of these values are 

Mco ^ - 1 

s “ jj, ! i)e]2/3 



(See appendixes D and E.) In transonic perturbation theory the terms 
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in 6 on the right-hand side are neglected. This can lead to fairly 
large eirors for even moderately large values of 6 since the approach 
to 0=0 is nonlinear and 


di 


“S 


dl. 


d0 


and 


PM 


d0 


.00 as 0 


(59) 


Judging from this one might expect, that quantitative agreement of tran- 
sonic pert\irhation analyses with experiment would not he so good. How- 
ever, in comparing two similar shapes with only sli^tly different 
thickness ratios hy transonic similarity considerations one would 
expect fairly good agreement. 


EXPERIMENTAL RESULTS 
Flow Field near 10° Wedge 


Figure 9 shows interferograms of the flow past the 10° semiangle 
wedge for l4 Mach numbers from O.7OO to O.892 and 1.207 to 1.465 (the 

1 ° 1 ° 

interferograms for the were very similar and hence are 

not shown here). Notice that the lines of constant density in the 
subsonic-flow interferograms are roughly elliptical in shape as predicted 
by the theory ( see appendixes A and B) . A supersonic flow region was 
first detected between Moo = O.7OO and 0.794- (the sonic line is shown 
as a dashed line in the figs.) and a shock emanating from the comer 
appears in the supersonic zone at Moo = 0.794. As the Mach number was 
increased, this zone grew larger and a shock appeared at the rear of it, 
while the shock emanating from the comer weakened and disappeared. 

This rearward shock was of the typical X type associated with a laminar 
boimdary layer, and the interferograms clearly indicated the separation 
of the boundary layer ahead of this shock. The similarity between the 
flow field at Moo = 0-892 and at M»o = 1-207 (figs. 9(d) and 9(e)) is 
striking; the base of the rearward shock has moved quite far back on the 
wake of the blunt trailing edge at Moo = 1.207 but in the vicinity of 
and ahead of the sonic line the two fields axe nearly identical except 

for the detached shock wave which appears about 1^ chord lengths ahead 

of the wedge at Moo = 1.207. As the Mach number was increased 
above 1.207, the detached shock moved in closer to the leading edge and 
finally "attached" at a Mach number quite close to the theoretical 
attachment Mach number of Moo= 1.4l8. Notice that the process of 
attachment is very continuous. The effect of the boundary layer is 
quite noticeable in the last few interferograms: This can be rou^ly 

accoxinted for by considering the boundary layer to change the shape of 
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the body by its displacement thickoess and then considering a nonviscous 
flow past this revised shape. On the wedge the boimdary layer will not 
grow so rapidly as on a flat plate because of the favorable pressure 
gradient and, in fact, the effect of the strong expansion around the 
corner is known to cause an almost complete collapse of the boixndary 
layer there. As the bow shock wave gets close to attachment the veloci- 
ties in the subsonic region behind it are getting very close to sonic 
velocity and hence the flow in this region is very sensitive to any 
sli^t curvature of the "revised shape" of the wedge. This accoimts 
for the shift of the base of the sonic line forward to the leading edge 
as the shock approaches attachment. The nonviscous theory would indi- 
cate that the sonic line would always begin at the comer and, at a 
Mach number just slightly above the shock- attachment Mach nmnber, the 
whole subsonic region would become sonicj then, with increasing Mach 
number, the flow behind the shock would be completely supersonic. As 
observed, the boundary- layer effect is to make the wedge have a curved 
surface and the sonic line actually moves slowly from the comer to the 
nose. Even with attached shock wave at M^ = 1.465 the flow behind 
the shock is not quite uniform (as nonviscous theory woiild indicate it 
should be) because of the effecti-f^e curved surface caused by the 
boxindary layer. 


Local Mach Nimober Distributions on Three Thin Wedges 

Figure 10 shows the variation of local Mach number distribution on 
1 ° 1 ° 

the surfaces of the ^ "^2 ^ semiangle wedges with free -stream 

Mach nijmber. This should be compared with figure 5 which shows the 
corresponding theoretical curves in terms of the transonic similarity 
parameters . The general behavior of the theoretical and experimental 
curves is quite definitely in good agreement. Particularly noteworthy 
is the slow variation of the local Mach number distribution near free- 
stream Mach number 1. 


Pressure-Coefficient Distributions on Three Thin Wedges 

The slow variation of the Mach number distribution in the range 
near M„ = 1 is obscured when the results are plotted in terms of pres- 
sure coefficient, since the pressiire coefficient changes a great deal if 
local Mach number is constant while the free- stream Mach number changes. 
A better parameter for presenting transonic pressure distributions would 
be p/pq (p/Pq' in case of a detached shock). Typical Cp distribu- 

1 ° 

tions are shown in figure 11 for the wedge (the results for the 10° 
and wedges were very similar and hence they are not presented) . 


4 
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The points shown were where the fringes intersected the body in the 
interf erograms . Since for a wedge the drag coefficient is proportional 
to the average Cp, the drag rise is evident in the subsonic distribu- 
tions as the point where Cp = 0 moves rearward with increasing free- 
stream Mach mmber. Linearized subsonic theory (which predicts Cp = 0) 
locates the Cp = 0 point at x/c = 50 percent. Figure 7 shows theo- 
retical reduced Cp distributions at various reduced free- stream Mach 
numbers. Again the qualitative agreement of these curves with experi- 
ment is evident. 


Shock-Detachment Distance for Three Thin Wedges 


Figure 12 shows the shock- detachment distance against reduced free- 
stream Mach number for the three thin wedges and includes the theoreti- 
cal values from reference 21. Here Vincent! and Wagoner's values for 


I 

-*00 


have been multiplied by 



1 ° 

0 = 7 ^ 

2 

1 ! 
“A 

0 

0 

11 

CD 


in order to make the transonic 


perturbation value of detachment reduced Mach number agree with the 

1° 

value from oblique- shock theory for the wedge: The reason for this 


was discussed in the section "On Comparing Theory and Experiment," 
namely, the difficulty of comparing transonic perturbation theory quan- 
titatively with experiment. Notice how rapidly the shock wave moves 
away from the wedge as the Mach number is decreased toward 1. 


Drag-Coefficient Variation with Mach Nimiber for Three Thin 

Wedges 

It was shown in reference 25 that the viscous effects on the wedge 
tend to compensate each other at the leading edge and the shoulder so 
that the over-all pressure drag is nearly the same as if the flow were 
inviscid. Thus it would be expected that the pressure-drag coefficients 
obtained by integrating the experimental pressure distributions would 
check the inviscid transonic perturbation theory. The reduced drag 


I 4 . *1 o 

In terms of Mach nmber, for the 7^ wedge the shock theory pre- 
dicts attachment at Mo^ = 1.33 = 1 * 68 ), while the transonic 

perturbation theory predicts M^o^ = 1.25 “ 1.19^- 
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coefficient used here was 



which is, in essence, the reduced drag coefficient of the upper (or 
lower) half wedge. This was done since the wedge model was regarded as 
the front half of a dbuhle-wedge profile and hence the value given here 
is the part of the reduced drag coefficient contributed by the front 
half of such an airfoil as in equation (42)], based on the chord 

of the double -wedge profile, which would be twice the chord of the model 
used here. Of course, this viewpoint is valid only for supersonic free- 
, stream Mach numbers. 

Figure 13 shows the reduced drag coefficients for the three thin 
wedges plotted against reduced Mach number. It is seen that the results 
give nearly a universal curve, which they should if the transonic simi- 
larity law is true, but that there are systematic variations with wedge 
angle. This is to be expected based on the discussion of the section 
"On Comparing Theory and Experiment." The vertical lines through the 
experimental points indicate estimated accuracy of the data. This fig- 
tire should be compared with figure 8, the theoretical reduced-drag- 
coefficient variation with reduced Mach number. It is obvious that the 
qualitative agreement of theory and experiment is good. In figure l4 
the theory and experiment are compared directly for the three thin 
wedges. Here the theoretical drag coefficients are shown with a verti- 
cal spread, the upper values for M^ > 1 corresponding to the use of 

-2(u - uj 

the pressure coefficient = rr and the lower values, to the 

P a* 

-2(u - uJ 

use of the pressure coefficient Cp = . The situation is 

vice versa for Moo <1* From this figure it is evident that the tran- 
sonic perturbation theory gives a good approximation to experiment. 


Flow Field at Moo = 1*44 for the 26.6° Wedge 

Figure 15 shows the experimental and theoretical constant- velocity 
lines in the subsonic region behind the detached shock wave for a 
26.57° semiangle wedge at Moo = 1.44. The theoretical analysis was 
made from relaxation calculations by Drebinger (reference 19) who used 
the flow equations with entropy variation behind the shock taken into 
account. The experimental constant-velocity lines were determined from 
the isopycnic lines of the interferogram by taking into account the 
lateral stagnation-pressure gradient behind the curved shock. The 
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isopycnic lines near this strong shock wave were probably slightly in 
error because of the "smearing out" of the pressure discontinuity across 
the shock in the side-wall boundary layers. It is seen that the agree- 
ment between theory and experiment on detachment distance and constant - 
velocity contours near the wedge is good. 

Figure 16 shows the surface pressure distribution from reference 17 
and the present experiments. Again it is seen that the agreement is 
good. 


Flow Field near the 8 . 8 -Percent Circular-Arc Section 

Figure 17 shows interferograms of the flow past the 8 . 8 -percent 
circular-arc section for l4 Mach numbers from O .718 to O .936 and 1.11 

to 1 . 500 . 

Supersonic velocity first occurred at = 0.825 (see the section 
"Critical Mach number") and in figures 17(c) and 17(d) a nearly symmetric 
supersonic zone is shown at M^o = 0,8^. No shock waves were apparent 
in this zone although a sensitive schlieren apparatus might have shown 
some weak shocks there. At M^ = O. 89 O the supersonic zone has grown 
rapidly and now terminates in the X- shock configuration. Further 
increase of the Mach number to M^„ = 0.935 (figs. 17(e) and 17(f)) 
shows the supersonic zone increasing laterally and the terminating shock 
moving rearward into the wake of the body. Figures 17(e) and 17(f) also 
show the density distribution at M^^ ~ 1.11 (the detached shock wave 
was just out of the field of view of the interferometer) ajid it is 
interesting to note the similarity between the flow field at Mj„ = 0.935 
and Mjjg Rj 1.11. It would appear as though the shock terminating the 
supersonic zone at M^o = 0.935 had moved rearward to form the trailing- 
edge shock (which is actually in the wake here because of the bl\int 
trailing edge) and the supersonic zone had grown laterally vmtil the 
sonic line joined with the detached shock far away from the body at 
M„ = 1 , thus causing an embedded subsonic zone in the supersonic flow 
with further increase in Mach number. 

With further increase of Mach number above Moo = 1.11, figures 17(g) 
to 17(n) show that the detached shock again approached the leading edge 
and the embedded subsonic zone decreased in size until finally the 
shock "attached" somewhere between Moo = 1.^0 and Moo = IA 50 (the 
theoretical value being Moo^ = 1.423). 


Local Mach Number Distributions on 8 . 8 -Percent Circular -Arc Section 

Figure I 8 shows the local Mach number distributions for the 
8 . 8 -percent circular-arc section as obtained from the experiments at 
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vsurious free- stream Mach numbers. Again it is apparent that the varia- 
tion of local Mach nmber distribution near M^ =1 is very slow and, 
indeed, the distribution for M^ = 1 could be interpolated from this 
figure with good accuracy. 

Figure 19 is a cross plot of the data of figure I8 except that here 
the data are given in transonic similarity parameters. This figure 
shows contours of constant reduced local Mach number on a plot of 
reduced free- stream Mach number against chordwise position. The dashed 
lines represent subsonic local Mach niimbersj the solid lines, supersonic 
local Mach numbers. Note again the slow variation of local Mach number 
distribution with free-stream Mach number near sonic velocity. 


Pressure-Coefficient Distributions on 8.8-Percent 


Circ\xlar-Arc Section 

Figure 20 shows the pres sure -coefficient distributions on the 
8.8-percent circular-arc section for various free-stream Mach ninnbers. 
The points shown are where the fringes intersected the body in the 
interferogram. Again the presentation in this manner obscures the 
interesting fact observed in figure 18. 


Drag-Coefficient Variation with Mach Ntimber for 8.8-Percent 

Circular-Arc Section 


Figure 21 shows the experimental determination of the drag coeffi- 
cient of the front part of the 8.8-percent circular-arc section. This 
again is of the nature of a fore-drag coefficient and, as shown 
in equation ( 43) , it should have a positive slope equal to 
^ "t 2 1 

c ’ y ~" i .^Dlw **00 = ■this is how the subsonic data have 


M»=l 


been Joined with the supersonic data. The vertical lines through the 
experimental points again indicate estimated acctiracy of the data.' For 
the case of an attached shock the pressvire distribution can be calculated 
using characteristics theory and the shock polarj however, a close 
approximation is obtained by considering the flow behind the shock wave 
to be Prandtl-Meyer flow. (This yields, approximately, parabolic- shaped 
bow and trailing-edge shock waves; see reference 33*) From this pressvire 
distribution the drag was calculated and is shown in figure 21. Taking 
into account the "reflected" characteristics from the shock wave would 
give more compression and increase the drag coefficient so that it would 
agree better with the experimental values of M^ = 1.4-50 and I.500 
shown in figure 21. 
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Note that the tests were made at low enough supersonic speeds to 
get definitely helow the drag-coefficient maximum at ;«1.20. 


Local Mach Number Distributions on a 12-Percent 
Biconvex Circular-Arc Airfoil 

• Figure 22 shows local Mach number distributions from reference 8 
for high-subsonic-speed flow over a 12-percent biconvex circular-arc ' 
airfoil (with turbulent boundary layer) . The data for the 8.8-percent 
circular-arc section at two supersonic speeds have been scaled according 
to the transonic similarity laws to the 12-percent case and are shown 
for the front half of the 12-percent airfoil in. figure 22. The back 
half for these two cases has been faired in using a Prandtl-Meyer expan- 
sion which should be approximately correct (a more accurate determina- 
tion could have been made using characteristics theory and the shock 
polar) . At M^ = 1.58, the theory indicates that the shock is attached 
with sonic speed just behind the shock on the leading ed^e, so that the 
distribution can be obtained by standard methods mentioned above; again 
the Prandtl-Meyer expansion approximation was used for the distribution 
at M^ =1.58 in figure 22. 

The behavior of the Mach number distributions is similar to those 
shown previously, except in this case the movement of the shock termi- 
nating the local supersonic zone is shown. Apparently little change in 

local Mach number distribution occurs between M = 0.936 and M =1.29 

00 00 ^ 


Drag-Coefficient Variation with Mach Number for a 
12-Percent Biconvex Circular -Arc Airfoil 

The data of figure 23 were converted to pressures which were inte- 
grated to give the pressure-drag coefficient for the various free-stream 
Mach numbers. In addition the drags of the front and back halves are 
shown separately. The drag- coefficient variation between Mo<, = O.96 
and Moo = 1.20 was based on constant local Mach number distribution at 
values interpolated between the curves for Moo = O.936 and = 1.29. 
The data were faired into the curves for attached shock wave calculated 
on the Prandtl-Meyer expansion basis. It is seen that the fore-drag 
coefficient has a maximum after Moo = 1 while the drag coefficient of 
the rear part has a maximum before Moo = 1* The over-all airfoil has a 
maximum drag coefficient just before Moo = 1 in order for the curve to 
have the slightly negative slope at Moo = 1 given by equation ( 4l) . 
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CONCLUSIONS 


✓ ✓ 

The transonic similarity theory of Von Karman and Guderley was 
checked and found to he in good agreement with experiment for thin 
wedge profiles near a free-stream Mach number of 1. . 

The results of theoretical calculations, using transonic perturba- 
tion theory, made by Guderley and Yoshihara, Vincenti and Wagoner, and 
Cole for a wedge in transonic flow were checked experimentally at high- 
subsonic and low- supersonic speeds for three wedges of different angles 
and were found to be in good agreement with experiment. 

The flow field and the surface pressvire distribution for a 26 . 6 ° 
semiangle wedge at a free-stream Mach number of l.i)-4 were obtained 
experimentally and were found to be in excellent agreement with the 
theoretical calculations of this flow made by Drebinger. 

The pressure distributions and drag coefficients for an 8 . 8 -percent 
circular-arc section followed by a straight section and for a 12 -percent 
biconvex circular-arc airfoil were presented completely through the 
transonic range. It was shown that some difficulty arises in comparing 
two-dimensional transonic perturbation theory with experiment, since 

this theory neglects thickness-ratio terms of order (t/c)^^^ and 
higher; for even moderate thickness -ratios this -will cause noticeable 
deviations from more exact theory. , , . 

It was shown from some physical arguments that the local Mach 
number distribution on bodies traveling through an infinite fluid has 
a stationary value at free-stream Mach mamber 1. This was verified 
experimentally for the case of two-dimensional flow. It was shown that 
this concept implies a drag- coefficient maximum just below free-stream 
Mach number 1 for all bodies in steady flight. This fact can be used 
to obtain the variation of local Mach number distribution on bodies 
completely through the transonic range of velocities from wind-tunnel 
tests, provided small models are used so that tests can be carried well 
above critical Mach number and to low-enough supersonic Mach numbers 
so that the bow shock wave is detached a chord length or so. 


California Institute of Technology 
Pasadena, Calif., June 1, 1951 
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APPENDIX A 

ASYMPTOTIC REPRESENTATION OF COLE'S SOLUTION FOR LARGE 
NEGATIVE VALUES OF REDUCED FREE-STREAM MACH NUMBER 


Cole's solution for the high-subsonic-velocity flow past a thin 
wedge (reference 22) is given as follows (in Cole's notation) 

l2zz-^^^ O" sinh X (v^ - v) 

y(z,vjzi) = dX 


( 61 ) 


/ \ l/^ ?/^ r~ cosh X ( v„ - v) 

x(z,v;zj . 1 - ' z ' 1 J 2 / 3 (Xz)J_j^j(Xzi)X dX 


(62) 

where the center line of the we 4 ge is at y = 0 ; the leading edge, 
at X = 0 ; the shoulder, at x = 1 ; and 




^3/2 2 

y V U 

) ■ 3 

-(7 + 1 ) 


3/2 




M. 




3/2 


u 


.3/2 


-(r + 1) 

a* 


Vq = (7 + i)e 


( 63 ) 


and the other notation is the same as that in the present paper. 

Using the standard methods of partial- fraction expansion one may 

write 


sinh X 
sinh 


. 1 . X + 2 f sin 


,2 2 

^ ^o 


X^v^2 ^ ^2^2 


( 64 ) 
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Equation (67) is from reference 3 ^} equation (66) is obtained by differentiation of 
equation (67) with respect to 7. 
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Figures 5 ^ 6 , and 7 were calculated from this equation for v = 0 , for the cases 
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1^3 

one can write equations (68) and (69) for large values of z and zj, 
as 




Eliminate v tetween equations (7^) and (75) for > z — 



and for z > Z]^ — simply replace x by x - 1 and z - Z]^ 
by zi - z in equation (76) . Thus the lines of constant Mach number 
are ellipses with center on y = 0, with ratio of semiaxes equal to 
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which is precisely the solution given by the linearized subsonic 
( see appendix B) . 

Now, in the notation of the present paper. 


- ^1) ^ 2n 
Vq 3 


(.1)3/2 . (.,j3/2 


and since I - loo is small on the wedge and since 

Cp = -2U-IJ 

one can write 




/^_ _ 
1 1 - 

3/2 

- 1 - 

L 

r 


Hence, 




n(z - z-^ 




1 3 '-p . 


- 1 


so, approximately. 


- zp) ^ \i/2 n Jt V 

V r'V ~~2 


1 - 


Q » ^*^1 


Z • Zn > 00 


Similarly, for large values of z and z^ it follows that 

.1/6 




\l/b 

I ~ Ij z. 


zi — >co for z - zp small 


theory 


(78) 


( 79 ) 


Substituting equations (78) and ( 79 ) into equation (16) one gets the 
exact linearized subsonic solution for constant-velocity lines (see 
appendix B) . Therefore on the wedge (y = O) , from equations (76), (78), 
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and (79), one has approximately for large values of z and 



Solving this for Cp 



or . 


C 


P 


-26 


n \fl 



loge 



(80) 


which is precisely the linearized subsonic solution for flow past a 
wedge (see appendix B) • Thus Cole's solution far away from = 1 
tends exactly to the linearized subsonic solution. 


/ 


/ 
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APPENDIX B 


LINEARIZED SUBSONIC AND SUPERSONIC FLOW PAST WEDGE 
AND CIRCULAR-ARC SECTIONS 
Linearized Subsonic Flow past a Wedge 


Let the wedge center line be on y = 0, with the leading edge at 
X = 0 arid the shoulder at x = Ck Then the Inconqpressible-flow problem 
is to find an analytic function u - iv such that v = 0 on y = 0 
except for 0 < x <.c where v = U6 and u - iv = 0 at infinity. 

Such a function is 


u 


IV = 

(z/c) - 1 


( 81 ) 


where 

z = X + iy 

U free -stream velocity 

Thus on y = 0, 



-2u 

U 



jt 


x/c 

1 - (x/c) 


( 82 ) 


Using the Prandtl-Glauert transformation, for linearized subsonic flow 


^1 - Mop 1 ~ (jt/c) 

or in transonic similarity notation 


(83) 





-2 


•«\/^ 


logp 


^ c 

1 - (x/c) 


( 84 ) 
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For the incompressible case the lines of constant pressiire in the fluid 
will be where 


z/c 


(z/c) - 1 

but these are circles with centers at 

e d 


- = Constant 


X = 


“«Cp ^ 


y = 0 


( 85 ) 


and radii 


1 - e 0 


Cp 

T 

e 

— nCp 

1 - e 


In the Prandt l-Glauer t transformation, the y distance is stretched by 

the factor \j 1 - Moo^ as is the pressure coefficient so the lines of 
constant pre ssure ( and hence density) are ellipses with ratio of axes 

equai to \/l - given by the equation 

r 




2 \/l - Moo^ sinh 4- — 




Mo 


2 

+ 


itC, 


\/r: 


M„ 


20 


X + 


^c^ \l 

1 -m2 

--P V 

00 


20 


„ Cp \/l\. 


sinh ii 
2 


= 1 (86) 


Linearized Supersonic Flow past a Wedge 

From the Ackeret theory the pressure coefficient in supersonic 
flow is proportional to the slope and for the wedge yields simply 


s = 


20 


- 


(87) 
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or 


2 


(88) 


Linearized Sixbsonic Flow past a Circular -Arc Section 

For the circular -arc section, the slope of the surface varies 
almost linearly with distance from the zero-slope point along the axis 
of the profile. For the section shown in figure 1 then, with the center 
line on y = 0, the leading edge at x = 0, and the zero-slope point 
at X = c, the Incompresslhle-flow problem is again to find an analytic 
function u - iv such that on y = 0, v = 0 except for 0< x< c 

where v = 2U — (l-§) where t is the half thickness at x = c and 
u - iv = 0 at infinity.. Such a function is 




log. 


l/c 


(z/c) - 1 


- II 


On the wedge (y = 0, 0< x< c), then 


r - 2si = it i 

Po " U " Jt c 




( 89 ) 


(90) 


so the linearized subsonic solution is 


Cp 


_ i i 

It c 


v/TT 


U 


(1 - |)loge 


A 


x/c 


1 - (x/c) 


+ 1 




(91) 


or 


Ct) = 


n\pTc 




(x/c) 


+ 1 


(92) 


The minimum Cp is obtained by differentiation, and one finds that 


dC 


d(x/c) 


2 -= 0 
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at the point where log- r- = -=- 

®l-(x/c) x/c 

of this transcendental equation is 


and numerically the solution 


I = 0.783 


which gives 


^Pmin “ 


-I.626 (t/c) 

\A - 


( 93 ) 


( 94 ) 


Linearized Supersonic Flow past a Circular -Arc Section 


The result here is again sln5)le from the Ackeret theory; 



( 95 ) 


(96) 



50 


NACA TN 2560 


APPENDIX C 

TRANSONIC SHOCK POLAR 


The equation of the shock polar in the hodograph plane is 


2 o* U 

^2 = (u - n) (97) 

^ + 2 IL ^ 

U 7 + 1 a* " a* 


where U is the velocity ahead of the shock and u and v are velocity 
con^jonents behind the shock parallel and perpendicular to the direction 
of U, respectively. Making the transonic approximation in this equa- 
tion let 


— * 
u = a + u 

V = V > 

U — a * 1 “ Uoo 


( 98 ) 


Substituting into equation (97) ^ neglecting higher powers of the per- 
turbation velocities, one obtains 


Letting 


v2 . 7+1 
2 


(Uoo - U)^ 



(99) 


u' = (7 + 1) 

a* 

v’ = (7 + 1) ^ 
a* 


( 100 ) 


one then has 


2(v’)2 = (u„’ - u')2 {u’ + uj) 


( 101 ) 
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The wedge arigle for detachment of the shock will now be 
given by the maximum value of v'. This is easily seen to occur 

at u* = - i Uoo’ giving v^^y ■ = i \/3 u.3/2. Since within the tran- 
sonic approximation 


on the wedge and 


v' = (7 + 1)0 


t 

U.00 




1 


( 102 ) 


this implies that 


(7 + l)0max ^ ~ jy2 (^“^ “ 1)^^^ (IO3) 

or, viewed in another li^t, this ingslies that the reduced attachment 
Mach number is . - 



1 .3 

[(7 + l)0]2/3 42/3 


1.19 


(104) 


Similarly, the wedge angle for obtaining exactly sonic velocity behind 
the shock is given by the value of v' where u* = 0. This is 

v' = -;== (uoo*)^/^. Again using equation (102) this can be written: 


(7 + 1)00 sj-^(Meo^ -1)3/2 

/2 


(105) 


and, viewed in another light, this implies that the reduced Mach nmber 
for which sonic velocity is obtained behind the shock on a wedge is 


^cog = - V., = 2V3 ^ 1.26 


(7 + 1)0 ^/3 


(106) 
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APPENDIX D 

VARIATION OF REDUCED MACH. HUMBER AT WHICH SONIC VELOCITY OCCURS 
BEHIND AN OBLIQUE SHOCK WITH FLOW -DEFLECTION ANGLE 


The ohllque-shock relations can he vritten 


M2^ = 


7-1 2 

1 + L— Mi^ 


M^2 cos^P 

7 Mi 2 sin^p - / - “ ■ 1 + Mi^ sin p 


where 


Mn 


M2 

P 

6 


tan 6=2 cot P 


M]_^ sln^p - 1 
Mi^( 7 + cos 2P) + 2 


Mach number ahead of shock 
Mach number behind shock 
shock-wave angle 
flow deflection 


For M 2 = 1, eliminating p between these two equations yields 


tan 6 = 


f(Ml) 


Ii^ - 1 + f (Ml) 


27 


TM]^^ - f(Mi) \l . 7+1 


1 + 


27 


f(Mi) 


where 


f (Ml) = _ 1 + j|l + (Mi2 - 1) + 


(107) 

(108) 


(109) 
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Expanding the ri^t-hand side in powers of - 1 (assumed small), 

one obtains 


(7 + 1) tein 6 = 


(Mjg , 1)3/^ 

gl/2 


^ " T ' <r - 1) + . . . 

16(7 + 1) 


(no) 


Reverting this series and letting 0 ~ tan 0 and M]l = Moog> one finds 


. 1 


’s [(r + i)e]2/3 


. pl /3 


1 + .. , 197 -t.g, e2/3 + o(eV3) 
',1/3 

2 




(111) 


For 7 = l.i)-. 


197 + 2 

WJ- 


73 


= 1.122 


i 
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APPENDIX E 

VARIATION OF REDUCED MACH NUMBER BEHIND A PRANDTL-MEYER 
EXPANSION FROM M = 1 THROUGH AN ANGLE 6, WITH 6 

The exact relation here is 




PM 




.-1 i/ij 

TM 


Expanding the right-hand side in terms of 




- 1, using 


(112) 


-1 ^^' 4 . 

tan X=X--:^+-^--;:r+* • • 
3 5 7 


one obtains' 


(7 + 1)0 = |( hp „‘ 


2 , V 


3/2f-(-l)\(MpM^ - 




-f 


( 113 ) 


where 


7 + 1 

^n = — F" 


( 7 - 1 
■ V7 + 1 


^n+l 


Reverting this series the first few terms are 

3/2 


^M 


- 1 = 


3(7 + 1) 
. 2 


7 5(7*'+ 1) 


3(1^ , o(a‘‘/3) 


J 


(ll4) 
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Therefore, 
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Flgvire 2.- Typical finite-fringe interferogram. 8.8-percent circular- 

arc section at M = 1.200. 
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Figure 3*- Typical superimposed finite-fringe interferogram. 10^ 

angle wedge at M = 1.278. 


semi- 



Figure 4.- Increment in Mach number per fringe contour against local Mach 
number for air (7 = 1.4), For these tests X/kZp = 0.025 + 0.001. 





Figure 5-- Theoretical reduced local Mach number distributions on a 
wedge near Mach number 1. Data from references 20 to 22 and 29. 











Figure 6.- Theoretical reduced local Mach number against reduced free- 
stream Mach number at several chordvise stations on a wedge. Data 
from references 20 to 22 and 29* 
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Figlire 7.- Theoretical reduced pressure-coefficient distributions on a 
wedge near Mach number 1. Data from references 20 to 22 and 29. 
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(a) 10° semiangle wedge. 

I 

Figure 10.- Local Mach number against x/c for increasing free-stream 

Mach number. 
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Figure 10.- Continued. 
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Figure 10 .- Concluded. 
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Figxire 12 .- Shock-detachment distance against reduced Mach number for a 

wedge . 
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Figure 13 •- Experimental reduced drag coefficient against reduced Mach 

nvHDiber on a wedge. 
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Figure 17- - Continued. 





Figure 1T>- Continued 
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Figure I8.- Local Mach number against x/c for increasing free-stream 
Mach manber. 8 . 8 -percent circular-arc section. 
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SIMILARITY EXTENSION OF PRESENT TESTS 
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Figure 23.- Drag coefficient against Mach number for 12-percent biconvex 
circular-arc airfoil. (From data of fig. 22 on basis of constant 
Mach number distribution from M^ = 0.96 to 1.20.1 












